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Abstract 

The radiative capture process n+p — ► d+j provides clear evidence for meson 
exchange currents in nuclear physics. We compute this process at low energies 
using a recently developed power counting for the effective field theory that 
describes nucleon-nucleon interactions. The leading order contribution to this 
process comes from the photon coupling to the nucleon magnetic moments. At 
subleading order there are other contributions. Among these are graphs where 
the photon couples directly to pions, i.e. meson exchange currents. These 
diagrams are divergent and require the presence of a local four-nucleon-one- 
photon counterterm. The coefficient of this operator is determined by the 
measured cross section, <r expt = 334.2 ± 0.5 mb, for incident neutrons with 
speed |v| = 2200 m/s. 
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I. INTRODUCTION 



The radiative capture n+p^d + 'yisa classic nuclear physics process where meson 
exchange currents play a role. For protons at rest and incident neutrons, with speed |v| = 
2200 m/s, the cross section for this process has the experimental value, <j expt = 334.2±0.5 mb 
Naively, one expects that an effective range calculation of this cross section |||| would 
be very close to this. However, such a calculation gives a value which is approximately 10% 
smaller than <j cxpt . As first suggested by Brown and Riska PJ, this discrepancy can at least 
partly be accounted for by the inclusion of meson exchange currents. More recent work 
by Park, Min and Rho || using effective field theory with Weinberg's power counting |5| 
for the nucleon-nucleon interaction, a resonance saturation hypothesis for the coefficients of 
some operators and a momentum cutoff finds the value a = 334 ± 2 mb. This prediction is 
relatively insensitive to the value of the cut-off and is compatible with cr expt . 

Recently, a consistent power counting for the nucleon-nucleon interaction has been es- 
tablished |7]]. At leading order in Weinberg's scheme, pion exchange is included in the NN 
potential and it is interated to all orders to predict the NN scattering amplitude f8|. How- 
ever, the work of [|7|,PJT0| shows that iterating the pions without including the effects of 
operators with explicit factors of the quark masses or derivatives does not give a systematic 
improvement in the prediction for the NN scattering amplitude. Using the power counting 
of 0, the bubble chain formed by multiple insertions of the momentum independent four- 
nucleon operator gives the leading scattering amplitude for systems with large scattering 
lengths. Higher derivative operators, operators involving insertions of the light quark mass 
matrix and pion exchange are of subleading order and are treated in perturbation theory. 
The expansion parameter used in |7|] is Q ~ |p|,m,r, and one expands in Q/A while keep- 
ing all orders in aQ (see also [|ll|]). Here A is a nonperturbative hadronic scale and a is 
the scattering length. At next-to-leading order (NLO) in the Q expansion simple analytic 
expressions can be derived for physical quantities. Various observables in the two-nucleon 
sector have been determined at NLO with this power counting, such as the electromagnetic 



form factors and moments of the deuteron [12], the polar izabilities of the deuteron [ 13 



Compton scattering cross sections ] T4| , |To' ] and parity violating observables pB|Ji7[ ]. 

In this work we compute the cross section for the radiative capture of extremely low 
momentum neutrons n + p — > d + 7 at NLO in the effective field theory Q expansion. 
Capture from the 3 Si channel is suppressed in the expansion compared to capture from the 

channel. At zero-recoil, the amplitude for capture from the 3 S\ channel vanishes since 
it is simply the overlap of two orthogonal eigenstates of the strong Hamiltonian. Hence, at 
leading order, only the isovector 1 So capture occurs and it arises from the isovector magnetic 
moment interactions of the nucleons. Since the amplitudes for capture from the and 3 Si 
channels do not interfere, at NLO the cross section comes only from the amplitude for 
capture from the ^So channel. At this order there are contributions from a single insertion of 
four-nucleon operators with two derivatives, from a single insertion of four-nucleon operators 
with an insertion of the quark mass matrix, and from the exchange of a potential pion. The 
potential pion exchange occurs in graphs with a nucleon magnetic moment interaction and 
also in graphs where the potential pion is minimally coupled to the electromagnetic field. 
The later contributions are historically called meson exchange currents. In addition to these 
contributions there is also a contribution from a four-nucleon-one-photon contact interaction. 
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The coefficient of this operator has not been previously determined and a major purpose of 
this paper is to fix its value. 



II. EFFECTIVE FIELD THEORY FOR NUCLEON-NUCLEON INTERACTIONS 



The terms in the effective Lagrange density describing the interactions between nucleons, 
pions, and photons can be classified by the number of nucleon fields that appear. It is 
convenient to write 

C = C + d + C 2 + ..., (1) 

where C n contains n-body nucleon operators. 

Cq is constructed from the photon field A^ 1 = (A , A) and the pion fields which are 
incorporated into an SU{2) matrix, 



2!in T1 _(*°/s/2 



E = exp m , n^-r_; /v ,j , (2) 

where / = 131 MeV is the pion decay constant. £ transforms under the global SU(2) L x 
SU(2) R chiral and U(l) em gauge symmetries as 

E -> LT,R\ £ -> e iaQem J:e' iaQem , (3) 

where L e SU(2) L) R e SU(2) R and Q e m is the charge matrix, 

«-=(JS) • < 4 > 

The part of the Lagrange density without nucleon fields is 

C = i(E 2 -B 2 ) + ^TrD^D^ + ^ATrm^E + E 1 ') + ... . (5) 
2 o 4 

The ellipsis denote operators with more covariant derivatives D^, insertions of the quark mass 
matrix m q = diag(m u , m^), or factors of the electric and magnetic fields. The parameter A 
has dimensions of mass and m 2 = \{m u + m d ) = (137 MeV) 2 . Acting on E, the covariant 
derivative is 



= <9 M £ + ie[Q em , T]A^ . (6) 

When describing pion-nucleon interactions, it is convenient to introduce the field £ = 
exp (in//) = VS. Under SU{2) L x SU{2) R it transforms as 

i - L£rf = USR\ (7) 

where U is a complicated nonlinear function of L, R, and the pion fields. Since U depends on 
the pion fields it has spacetime dependence. The nucleon fields are introduced in a doublet 
of spin 1/2 fields 
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N 



(8) 



that transforms under the chiral SU(2)l x SU(2)r symmetry as N — > UN and under the 
U(l) em gauge transformation as iV — > e ta ® em N. Acting on nucleon fields, the covariant 
derivative is 



D^N = (0„ + V I1 + ieQ em A^)N , (9) 



where 



V^^D^ + gDpt) = ^(Zd^ + gdrZ + ieA^QU-ZQem?)) ■ 

The covariant derivative of N transforms in the same way as N under SU(2)l x SU(2)r 
transformations (i.e. D^N — > UD^N) and under U(l) gauge transformations (i.e. D^N — > 
e iaQ * m D^N). 

The one-body terms in the Lagrange density are 



d = rt \iD Q + 



2M 



+ WH^ ( k + ^[£V 3 £ + £t 3 ?] )tr-BN + ... 



2M 



where M = 939 MeV is the nucleon mass and kq 



(10) 



2 (^%> ^V. 



0.4399 and «i 



o(/«p — K n ) = 2.35294 are the isoscalar and isovector nucleon magnetic moments in nuclear 
magnetons. The nucleon matrix element of the axial current is qa = 1.25. 
The two-body Lagrange density needed for NLO calculations is 



+ 



Cf 5l) + £>f Sl) ATrm g ) (N T (N T P t N) 
(iV T P i JV) t (n t \Pif) 2 + & 2 Pj - 2£)PE> 



N) +h.c. 



- (c£ So) + D£ So) \Trm q ) (^P^^^PiN) 



+ 



a 



(N T PiNY [N 



JT 



Pi 3 +15 Pi- 2 i DP i T) 



N ) +h.c. 



+ 



eh x (N T P i N)\N T P z N)B l - eL 2 %e ijk (N T P,N)\N T P j N)B k + h.c. 



where p and P« are spin-isospin projectors for the spin-triplet channel and the spin-singlet 
channel respectively, 



P 

P; 



Tr P/P, 



Tr PIP,- 



1 
1 



(12) 



The cr matrices act on the nucleon spin indices, while the r matrices act on isospin indices. 
The local operators responsible for S — D mixing do not contribute at NLO. Terms in £2 
involving the pion field have been neglected in eq. ([TT|) . 
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The values of the coefficients of the four-nucleon operators in £2 depend on the regular- 
ization and subtraction scheme that is adopted. The power counting of 0] is manifest in the 
power divergence subtraction scheme, PDS, and we shall use it in this paper. (Momentum 



subtraction schemes can also have this power counting jFJIyJ). In PDS one works in D 
dimensions and the poles at both D = 4 and D = 3 in the loop integrations associated 
with Feynman diagrams are subtracted. The D = 4 poles are from logarithmic ultraviolet 
divergences and the D = 3 poles are from linear ultraviolet divergences. There is some free- 
dom in how the Lagrangian is continued to D-dimensions. We choose to keep the Pauli spin 
matrices three dimensional and continue the derivatives to D-dimensions. This is similar 
to the scheme proposed by t'Hooft and Veltman |20| for chiral gauge theories and ends up 



being convenient since the n + p — > d + 7 amplitude is proportional to the antisymmetric 
Levi-Civita tensor, e^. At NLO two basic divergent integrals are encountered. The first is 

^y-D , d P-i) q ( 1 \ 



2) J (2tt)( d - 1 ) \q 2 + a 2 J 

" [Va ] L { 2 ) (4tt)(°- 1 )/ 2 ' {6) 
Iq has no pole at D = 4 but does have a pole at D — 3. Its value in the PDS scheme is, 

/ o PDS =(£> ( "-^ ) ' (14) 



The second is the two loop integral, 

1= \2J J (2ix) D - 1 (2u) D - 1 q 2 + a 2 FT^ (q - l) 2 + c 2 ' ^ ' 

Ii has no pole at D = 3 but does have the pole, — l/327r 2 (D — 4), at D — 4. Therefore I\ 
has the same value in minimal subtraction (MS) as in PDS f2~H| , 



T S = It" ^ — (>o g ( —7— ) + (16) 

where 

5= 1(^-1 -log Q)) (17) 

and is Euler's constant. Note that because of its logarithmic divergence, [3Ii/(D — 
1)] PDS = l[ DS + 1/967T 2 . There is considerable freedom in the precise way the subtractions 
are handled. For example, if the poles in D = 4 are subtracted with MS then 5 = —1/2. 
Finally, we stress that one cannot blindly evaluate the integrals in D-dimensions and subtract 
the poles to get the required PDS value. For example, if a and b are set to zero then ii has 
a double pole at D = 3. However, this is associated with a logarithmic infrared divergence 
in three dimensions, not an ultraviolet divergence, and so it is not subtracted. 

Most of the coefficients in £ 2 have been determined. At NLO the deuteron magnetic 



moment [T2|] is 
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e / 2M 7 (/i- 7 ) 2 \ , . 

/id= 2M \ Kp + Kn + n J ' ^ 

where 7 = \[MB with B = 2.225 MeV the binding energy of the deuteron and fi is the 
subtraction point. The coefficient L2 depends on the subtraction point in such a way that 
the physical quantity fid is fi independent. The experimental value of the magnetic moment 
of the deuteron is fid = 0.85741 nuclear magnetons and comparing this with the prediction 
above implies that the coefficient L2 (renormalized at /i = m n ) is, 

L 2 (m 7T ) = -0.149 fm 4 . (19) 

Note that, 

^P^N = ie i3k N T P k N, (20) 



and so the operator with coefficient L2 in eq.flTTD is the same as in \T2 

The coefficients of the four-nucleon operators in eq. (|TT|) that don't involve the electro- 
magnetic field have been fixed from comparison with experimental data on NN scattering. 
We will review this in the following section. The only unknown coefficient that contributes 
at NLO is L\ and it will be determined in this work. 



III. S-WAVE NN SCATTERING 



The 1 5'o NN scattering amplitude, at center of mass momentum p, has the expansion 
A^ So \p) = ££L_i A£ So ^(p), where A^ So \p) is of order Q n . At leading order only the four 
nucleon operator with no derivatives need be included and 



-a 



l + c£ So) M(fi + ip)/4ir 



(21) 



It is convenient to break the next order contribution into several pieces, 



Mi) Ail) .{iv) 



and, using PDS, ref. M found, 



All I) 



^2 P 



9l 



Q-s ) 



All 



a 







^2/ 2 J \ Ap 2 

g\ ( m n MA ( ^ o) 
P 



mi , ( Ap 

\ m: 



Att 



(/i + ¥) , ^TT 

m n 2p 



tan 



g\_ ( m^MAll 
2f 2 



An 



fi + ip 



+ i tan 1 



2p 

m n 

2p 



In 1 + 



Ap 2 



1 \ m*j 
_ I i n ( m " + 4 P 2 



1 
6 



A^ 



A es 0) 



('So) 



(22) 
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Actually, the above expression for A IV ^ is slightly different from that in because there 

some terms that appear above were absorbed into a redefinition of D\ . 

The scattering length gets contributions from each order in the Q expansion. To use these 
results for the scattering amplitude over a region that includes very low p it is necessary 
that the the leading order amplitude give almost the correct the scattering length. This can 
be achieved at NLO by reordering the expansion in the following way (T2| . Write 

cf s °) = c£ So) + Ac£ So \ (23) 



treat Cq nonperturbatively and ACq So ^ as a perturbation. Then in eqs. fl2T|) and (22) 
the following changes occur, 

cf o) - C§*» , D^ml^D^ml + AC^ So \ (24) 

— /in \ 

The coefficient Cq is no longer independent of the light quark masses and is chosen so 
that the physical value of the scattering length is close to the first term in its Q expansion, 




,2 



+ 1 - 1 + Ji_ 1 (25) 
4vr/ 2 IT ' Md 1So) J K ' 2 "" \ n J 2 12 2m 2 ) ' K ' 



The subtraction point is arbitrary and the coefficients Cq So \ C^ So ^ and Sq W^+ACo 
depend on \x in such a way that A_° and A are independent of \x. These coefficients 
are determined from the measured NN phase shift. We will only need the first two of them 
and for these a fit over the region 7 MeV < p < 100 MeV finds, at fi = [1 



C ( 1So) {m n ) = -3.529 fm 2 , C^ S °\m^) = 3.04 fm 4 . (26) 



Discussions of the results of different fitting procedures can be found in |7|, p!2]jr8| , |22|j23|| 

In the 3 Si channel, identical formulae hold for the scattering amplitude and scattering 
length once the replacement 1 S'o —>- 3 Si is made for the superscripts. However, the fit to 
the data is done a little differently. For processes involving the deuteron it is convenient to 
constrain Cq So ^ so that A[_i^ gives the correct deuteron binding energy, B = 2.2255 MeV. 
This implies that 



where 7 = a/ MB. In this channel a constrained fit to the NN phase shift yields 

C { Q Sl \m % ) = -5.708 fm 2 , cf 5l) (m^) = 10.8 fm 4 . (28) 
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FIG. 1. Graphs contributing to the amplitude forn + p — > d + j at leading order in the effective 
field theory expansion. The solid lines denote nucleons and the wavy lines denote photons. The 
light solid circles correspond to the nucleon magnetic moment coupling to the electromagnetic field. 
The crossed circle represents an insertion of the deuteron interpolating field which is taken to have 
3 S*i quantum numbers. 



IV. CROSS SECTION FOR RADIATIVE CAPTURE 



The amplitude for the radiative capture of extremely low momentum neutrons n + p 
d + 7 has contributions from both the and 3 Si NN channels. It can be written as 



iA(np^d~f)=eX N t t 2 ct 2 [<r ■ k e(d)* ■ e( 7 )* - <x ■ e( 7 )* k-e(d)*]N 
+ ieY e ijk e{dy* k j e( 7 ) fc * (N T r 2 r 3 a 2 N) , 



(29) 



where e = |e| is the magnitude of the electron charge, N is the doublet of nucleon spinors, 
e( 7 ) is the polarization vector for the photon, e(d) is the polarization vector for the deuteron 
and k is the outgoing photon momentum. The term with coefficient X corresponds to 
capture from the 3 S\ channel while the term with coefficient Y corresponds to capture from 
the channel. For convenience, we define dimensionless variables X and Y, by 



mV 7 3 



Y = % ^ [^Y 

mV 7 3 



(30) 



Both X and Y have the Q expansions, X = X + X 1 + and Y — Y + Yj. + where a 
subscript n denotes a contribution of order Q n . The capture cross section for neutrons with 
speed |v| arising from eq. ( P9|) is 



a 



87ra 7 
M 5 lv| 



2\X\ 2 + \Y\ 



(31) 



where a is the fine-structure constant. 

At leading order, X and Y are calculated from the sum of Feynman diagrams in Fig. ([[]) 



and from wavefunction renormalization associated with the deuteron interpolating field |L2 



giving 



X 



«o 1 + 



7 M 

47T 



Y 



Ki 1 + 



47T 



(32) 



S 



where A[_f°\p) is the leading, order Q^ 1 , contribution nucleon-nucleon scattering amplitude 
in the ^So channel at an center of mass momentum p. The scattering length is related to the 
nucleon-nucleon scattering amplitude at zero momentum 

A<*>\0) = -^a™ , (33) 

and the experimental value for the 1 So scattering length is, a^ s °^ = —23.714 ± 0.013 fm. An 
analogous expression holds in the 3 Si channel. At leading order, A[_i°\o) = —Ana^^/M 
and A^\o) = -An/M-f. Using this in eq. fl32|) gives, 

X = , % = Ki (l - 7 a (1So) ) , (34) 

which implies the radiative capture cross section, 

ro 87ra 7 5 ^(a( l5 »))2 / 1 \ 2 

o LO = J - 1 — = 297.2 mb . 35 

|v|M 5 \ jd s °) J v ; 

This agrees with the results of Bethe and Longmire |2||| when terms in their expression 
involving the effective range are neglected. Eq. (^) is about 10% less than the experimental 
value, cr expt = 334.2 ± 0.5 mb [|]]. Because Xq vanishes we only need compute Y\ to obtain 
the cross section at NLO. 

At NLO there are contributions from insertions of the D 2 , Ci operators and from the 
exchange of potential pions, with the photon coupling both minimally to the pions and 
to the nucleons via their magnetic moment. In addition there is a contribution from the 
L\ four-nucleon-one-photon operator. These can be divided into two categories, those that 
build up the ^So scattering amplitude at NLO, Aq So \ and those that do not. Writing, 

y x = y(rcscatt) + y-(Ca) + y^B) + y-( ff ,E) + y(L 1 ) > we find the grap h S in FigS. (|), @, ® and 

(|) give contributions 

f(r escatt) =Ki 7M^5 ) (0) 

An 



Y (°a) 



7 2 
~2 



c ' So)+c ' Si) a^\o) + 2 y*y-^ (i + ^) ( o) 



8nf 2 \ m7T + 2j An V7 \ m*J + 2-f J J 

" 12^ + 4^ (0) l v 2(m. + 7 ) +ln l v „ J 

y(Li) _ r ^2 -A-f O) (0) f o 6 N 

r - 7 ( i So) ( 3 5l) ■ W 

The NLO ^So scattering amplitude Aq So ^ was given in the previous section. Notice that the 
meson exchange current contribution Y^E) depends upon the renormalization scale \i. The 
graphs contributing to this term have a pole at D = A and require a subtraction. This is the 
reason for the logarithmic /^-dependence. It is canceled by the /i-dependence of the constant 
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FIG. 2. Graphs contributing to the amplitude for n + p — > d + 7 a£ subleading order due to 
the insertion of the C2 and D2 operators. The solid lines denote nucleons and the wavy lines 
denote photons. The light solid circles correspond to the nucleon magnetic moment coupling of 
the photon. The solid square denotes either a C2 operator or a D2 operator. The crossed circle 
represents an insertion of the deuteron interpolating field . The last graph denotes the contribution 
from wavefunction renormalization. 



L\ so that Y\ is independent of fi. It is interesting to note that the contribution from the C2 
operators, Y^ C2 \ is not \x independent either. This explicit /i-dependence is also canceled 
by the /z-dependence of L\. The renormalization group equation for the subtraction point 
dependence of L\ is 



d_ 



o °0 



48tt 2 / 2 ' 



(37) 



Note that this is quite different from the renormalization group equation, 



d_ 



(CI 



0. 



(38) 



that L 2 satisfies. 

There is a NLO contribution that we have not explicitly included, a one-loop correction 
to the magnetic moments of the nucleons that is of order m n /(4Trf 2 ). However, by using 
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FIG. 3. Graphs contributing to the amplitude for n + p — > d + j at subleading order due to the 
exchange of a potential pion with the photon coupling to the magnetic moment of the nucleons. The 
solid lines denote nucleons and the wavy lines denote photons. The dashed line denotes a pion. 
The light solid circles correspond to the nucleon magnetic moment coupling of the photon. The 
crossed circle represents an insertion of the deuteron interpolating field. The last graph denotes the 
contribution from wavefunction renormalization. 



the value K\ = 2.35294, which follows from the measured values of the neutron and proton 
magnetic moments, in Y , this effect has been taken into account. Similarly, using the 
measured value for a*- 50 •* in eq. fl3~4|) includes the effects of y( rescatt ). Demanding that the 
NLO expression for Y give the measured cross section implies that, 

Y (Ca) + f (^> B ) + y(^ E ) + f( L i) = 0.92. (39) 

In eq. using n = m n , 5 given by eq. (JT7|) and ^f o) = -Ana^/M yields, F (C2) = 0.38, 
y(«-. B ) = -0.33, and Y^^ = 0.60. Note that for Y^ there is a significant cancellation 
between the two terms in the square brackets of eq. (PBJ) . About two thirds of the discrepancy 
between the measured cross section and the leading order expression is made up from the 
meson exchange current contribution. Most of the rest comes from F (Ll) . Eq. © implies 
that 

Li(m^) = 1.63 fm 4 . (40) 
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>o- 

XD--CX^ 

FIG. 4. Graphs contributing to the amplitude for n + p — > d + 7 subleading order due to 
meson exchange currents. The solid lines denote nucleons and the wavy lines denote photons. The 
dashed line denotes a pion. The dark solid circles correspond to minimal coupling of the photon. 
The crossed circle represents an insertion of the deuteron interpolating field. 

The value of L\ is quite sensitive to the precise way that the poles at D = 4 are handled. For 
example if MS is used (i.e. 5 = —1/2) then Y^ W ' E ^ = 0.37 which gives Li(m w ) = 3.03 fm 4 . 

With L\ determined, all the counterterms in the strong and electromagnetic sector that 
occur at next-to-leading order in the effective field theory Q expansion are known. It is 
interesting to see that in this framework there is nothing special about meson exchange 
currents. They are simply one of the several contributions at NLO, occuring along with the 
strong interaction corrections to diagrams where the photon couples to the nucleon magnetic 
moments. 

V. CONCLUDING REMARKS 

We have computed the cross section for the radiative capture process n + p — > d + 7. At 
leading order we recover the effective range theory result (when the terms involving ro are 
neglected) which is about 10% smaller than the measured cross section . At NLO there are 
contributions from perturbative insertions of the C2 operators, the D2 operators, potential 
pion exchanges and from a previously unconstrained four-nucleon-one-photon counterterm 
with coefficient L\. In order to reproduce the measured cross section, <r expt , we find that 
Li(m n ) = 1.65 fm 4 . In more traditional approaches, meson exchange currents are required 
to explain the value of cx expt . In effective field theory, the meson exchange current graphs 
are divergent and require regularization. As a result, their contribution to the cross section 
is not unique and depends upon the choice of regularization scheme. In addition, a local 
counterterm is required to absorb these divergences and its value is a priori unknown and 
scheme dependent. Having determined the value of L\ from the radiative capture cross 
section, other processes arising from electromagnetic interactions such as deuteron breakup 
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FIG. 5. Local counterterm contribution to the amplitude for n + p — ► d + 7 NLO. The 
solid lines denote nucleons and the wavy lines denote photons. The solid circle corresponds to an 
insertion of the L\ operator. The crossed circle represents an insertion of the deuteron interpolating 
field. 

e + d^e' + n + p can be computed at NLO. Work on this is in progress. 
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